GLOBAL REGULARITY CRITERION FOR THE 3D NAVIER STOKES EQUATIONS 
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^_ ■ Abstract. In this paper wc provide a sufficient condition, in terms of only one of the nine entries of the gradient 

tensor, i.e., the Jacobian matrix of the velocity vector field, for the global regularity of strong solutions to the 
three-dimensional Navier-Stokes equations in the whole space, as well as for the case of periodic boundary 
conditions. 
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1. Introduction 
The three-dimensional Navier-Stokes equations (NSE) of viscous incompressible fluid read: 

du 

— -^Au+(u-V)u + Vp = 0, (1) 

V-u = 0, (2) 

u(x 1 ,x 2 ,x 3 ,0) = u (x 1 ,x 2 ,x 3 ), (3) 

• ■ where u = (tti, u 2 , u 3 ), the velocity field, and p, the pressure, are the unknowns, and v > 0, the viscosity, is given. 

We set S/h = (d Xl , d X2 ) to be the horizontal gradient operator and A^ = d% x + d% 2 the horizontal Laplacian, while 
V and A are the usual gradient and the Laplacian operators, respectively. In this paper we consider finite energy 
solutions of the system (l)-(3) in the whole space R 3 , that decay at infinity. However, we remark that one can 
apply our proof, nearly line by line, to establish same result for the three-dimensional Navier-Stokes equations 
in a periodic domain. 

The question of global regularity for the 3-D Navier-Stokes equations is a major challenging problem in applied 
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analysis. Over the years there has been an intensive work by many authors attacking this problem (see, e.g., [6], 
[7], [9], [19], [21], [22], [24], [25], [32], [35], [36], [37] and references therein). It is well-known that the 2D Navier- 
Stokes equations have a unique weak and strong solutions which exist globally in time (cf., for example, [7], [19], 
[32], [35], [36]). In the 3D case, the weak solutions are known to exist globally in time. But, the uniqueness, 
regularity, and continuous dependence on initial data for weak solutions are still open problems. Furthermore, 
strong solutions in the 3-D case are known to exist for a short interval of time whose length depends on the 
physical data of the initial-boundary value problem. Moreover, this strong solution is known to be unique and 
depend continuously on the initial data (cf., for example, [7], [19], [32], [35]). 

Starting from the pioneer works of Prodi [28] and of Serrin [31], many articles were dedicated for providing 
sufficient conditions for the global regularity of the 3D Navier-Stokes equations (for details see, for example, the 
survey papers [21], [37] and references therein). Most recently, there has been some progress along these lines (see, 
for example, [2], [3], [10], [11], [13], [14], [16], [33], [34], and references therein) which states, roughly speaking, 
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that a strong solution u exists on the time interval [0,T] for as long as 

2 3 

ue L p (\0,T],L q ), with - + - = 1, for q > 3. (4) 

p q ' 

Moreover, that has also been some works dedicated to the study the global regularity of the 3D Navier-Stokes 
equations by providing some sufficient conditions on the pressure (cf. e.g., [3], [4], [5], [8], [17], [30], [39]). In 
addition, some other sufficient regularity conditions were established in terms of only one component of the 
velocity field of the 3D NSE on the whole space M 3 or under periodic boundary conditions (cf. e.g., [4], [15], [18], 
[26], [27], [38]). 

We denote by L q and H m the usual L q — Lebesgue and Sobolev spaces, respectively (cf. [1]), and by 

H0llg = ( / \4>( x )\ q dx 1 dx 2 dx 3 J , for every <f> G L q . (5) 

We set 

V = {<j> : the three-dimensional vector valued C£° functions and V • <j> = 0} , 

which will form the space of test functions. Let H and V be the closure spaces of V in L 2 under L 2 — topology, 
and in H 1 under H 1 — topology, respectively. Let uq G H, we say u is a Leray-Hopf weak solution to the system 
(l)-(3) on the interval [0,T] with initial value uq if u satisfies the following three conditoins: 

(1) ue C w ([0,T], H) n L 2 ([0,T],V), andd t ueL 1 {[0,T],V / ), where V is the dual space of V; 

(2) the weak formulation of the NSE: 



u(x, t) ■ 4>{x, t) dx — / u(x , to) ■ (f>(x , to) dx 
[u(x,t) ■ {4>t{x, t) + vA(fi(x, t))] dx ds 
[(«(#, t) ■ V)(/>(x, t)] ■ u(x, t) dx, 



for every test function <j> <E C°°([0,T], V), and for almost every t, to G [0, T]; 
(3) the energy inequality: 

\\u(t)\\l + " [ l|Vu(*)||!<fa<K*o)||!, (6) 

J to 

for every t and almost every to- 
Moreover, if uq G V, a weak solution is called strong solution of (l)-(3) on [0,T] if, in addition, it satisfies 

MGC*([0,T],y)ni 2 ([0,T],i7 2 ), and d t ue L 2 ([0,T],i7). 

In this case, one also has energy equality in (6) instead of inequality, and the equality holds for every to- 

In this paper, we provide sufficient conditions, in terms of only one of the nine components of the gradient of 
velocity field, i.e., the velocity Jacobian matrix, that guarantee the global regularity of the 3D NSE. Specifically, 
if uq G V, and if for some T > and some k,j, with 1 < k,j < 3, we have 

-^i- eL /3 ([0,T],L a (R 3 )); when k ^ j, and where a > 3, 1 < j3 < oo,and - + - < ^— , (7) 

oxk a p 2a 

or 

^€LP([0,T},L a (R 3 )); where a > 2, 1 < < oo, and - + % < 3( " + 2) , (8) 

oxj a p 4a 
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where u — (u\, m 2 , U3) is a weak solution with the initial datum ua on [0, T], then u is a strong solution of the 31? 
Navier-Stokes equations which exists on the interval [0,T]. Moreover, u is the only weak and strong solution on 
the interval [0,T] with the initial datum uq. In particular, if (7) or (8) holds for all T > 0, then there is a unique 
global (in time) strong solution for the 31? NSE with the initial datum uq. 

For convenience, we recall the following version of the three-dimensional Sobolev and Ladyzhenskaya inequal- 
ities in the whole space M 3 (see, e.g., [1], [7], [12], [20]). There exists a positive constant C r such that 






■||9x^|| 2 2r ||M|| 2 2P \\9 X3 Ml 



<a 



Msi 



(K3)' 



(9) 



for every ip E H 1 ^ 3 ) and every r E [2, 6]. Observe that in case of periodic boundary conditions, one would have 
instead of (9) the following inequality 



< CMW7 



r (IIMh + \\i>h)~ (HMIh + Hh)^ {\\d,M* + UhY 

3(r— 2) 



lH!(0)' 



(10) 



for every ip E H 1 ^) and every r E [2,6]. Here, f2 is the periodic box [0,L] 3 . We remark that one can apply 
inequality (10) instead of inequality (9), and the methods presented in this paper to establish the same results in 
the case of periodic boundary conditions. The details of the proof in the periodic case are omitted. 



2. The Main Result 

In this section we will prove our main result, which states that the strong solution to system (l)-(3) exists on 
the interval [0, T] provided the assumption (7) or (8) holds. 

Theorem 1. Let uq E V, and let u = (ui,i«2,W3) be a Leray-Hopf weak solution to 3D NSE, system (l)-(3), 
with the initial value uq. Let T > 0, and suppose that, for some k,j, with 1 < k, j < 3, u satisfies the condition 
(7) or (8), namely, 



duj{s) 



dx k 



ds < M; when k ^ j, and where a > 3, 1 < j3 < 00, and \- — < , (11) 

a p 2a 



duj{s) 



i).i 



ds < M; where a > 2, 1 < f3 < 00, and \- — < . 

a p Aa 



(12) 



for some M > 0. Then u is a strong solution of 3D NSE, system (l)-(3), on the interval [0, T]. Moreover, it is 
the only weak solution on [0,T] with the initial datum uo- 

Before we prove the main Theorem 1, we show the following two lemmas. 
Lemma 2. 



/ 



/ g dx\dx2dx- c 



<cuitw*Ml_\\nT ii^ 2 /n 2 /r n& 3 /ii 2 /r 110112, 



(13) 



where 2 < r < 3. 



Proof. Observe, first, that it is enough to prove the inequality for functions 4>,f,g E C^°(]R 3 ) and then passing 
to the limit using a density argument. 
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/ g dx\dx2dx' c 



< C [ max |0| ( f fdxi) ( f g 2 dx 
Jk 2 Xl \Jr J \JM 



1/2 



dx2dX3 



<C 



<C 



max \cf>\ I dx2dx3 

R 2 V Xl 



1/r 



R 2 \JWL 



f 2 dx 1 



dx2dxs 



i :! 



g dx\dx2dx: 



1/2 



\<j)\ r \d Xl <j)\ dxidx 2 dx 3 



l/r 



1 1/r 

3- 



1/r 



<CW 2 — 11^011^11/H^ \\d x J\\ x 2 /J ||9x 3 /||2 /7 IM 



f r ~ 2 dx2dx-^ 



1 1/r 



cfcci 



lffl|2 



Lemma 3. 



where 2 < r < 3. 



/ 5 dx\dx2dx- c 



<c\\4>\\F IIS^fjT ll/ll^ 2 ll^/H^ R 2 /|| 2 /r IMh, 






Proof. Here, again, it is enough to prove the inequality for functions </>, f,g G C ( 
f g dxidx2dx3 <C I max|0| / f 2 dx 3 I / g 2 (ia;. 



<C 



<C 



I max |0| I dx\dx2 



1/r 



/ <ix3 I dxidx^ 



\4>\ r 1 \d X3 4>\ dxidx 2 dx 3 



1 1/r 



1/r 



1/r 



<cM 2 — R.tfimi/iir" ii^/iir ii9x 2 /n 2 /7 iiffii 



/ r - 2 dx\dx2 



1 1/r 



dxa 



dxi(ix2 



g dx\dx2dx- c 



9 2 



1/2 



D 



(14) 



□ 



Proof of the Theorem 1. Without loss of generality, we will assume that j — 3 and k = 1 in (11) and (12), 
namely, 



du 3 (s) 



dx\ 



P 



ds < M; where a>3, 1</3<oo, and 



3 2 a + 3 

a f3 ~ la 



(15) 



du 3 (s) 



dx 3 



rfs < M; where a > 2, 1 < /3 < oo, and — + - < — - -. 

a p 4a 



(16) 



It is well-known that there exists a global in time Leray-Hopf weak solution to the 3D NSE, the system (l)-(3), 
in whole space M 3 if u a e H (see, e.g., [7], [21], [23], [25], [32], [36]). It is also well-known that there exists a 
unique strong solution for a short time interval if u$ G V. In addition, this strong solution is the only weak 
solution, with the initial datum uq, on the maximal interval of existence of the strong solution. 

Suppose that u is the strong solution with the initial value uq E V such that u G C([0, T*), V)f)L 2 ([0, T*), H 2 ), 
where [0, T*) is the maximal interval of existence of the unique strong solution. If T* > T then there is nothing to 
prove. If, on the other hand, T* <T our strategy is to show that the H 1 norm of this strong solution is bounded 
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uniformly in time over the interval [0, 7"*), provided condition (15) or (16) is valid. As a result the interval [0, T*) 
can not be a maximal interval of existence, and consequently T* > T. Which will conclude our proof. 

From now on we focus on the strong solution, u, on its maximal interval of existence [0, 7"*), where we assume 
that T* < T. As we have observed earlier the strong solution u will also be the only weak solution on the interval 
[0,7"*). Therefore, by the energy inequality (6), for Leray-Hopf weak solutions, we have (see, e.g., [7], [32], [35] 
or [36] for details) 

Il«(*)ll2 + " / \\Vu(s)\\Us < Ki, (17) 

Jo 

for alH > 0, where 

Ki = Huolli- (18) 

Next, let us show that the H 1 norm of the strong solution u is bounded on interval [0, T*). 

2.1. || Vhit|J2 estimates. First we obtain some estimates of the horizontal gradient. Taking the inner product of 
the equation (1) with — A/jU in L 2 , we obtain 

= / [(u ■ V)u] • A/jU dx\dx2dx 3 . 

Jr 3 

By integration by parts few times, and using the incompressibility condition (2), we get 



(u ■ V)m • AftU dxidx2dx3 



3 2 

E^ du k duj du k 
/ — jr^a — dx\dx2dx^ 



k,j=l 1=1 



dxj dxi dxi 



<9«i\ 3 (du 2 \ 3 fdui du 2 \ 


f 9«i^ i ( du 2 \ x du\ du 2 


dxi) \dx2J \dx\ dx 2 ) 


\dx2) \dx\) ' 8x2 dxi 



E 

k,l=l 



du k du 3 du k 
dx 3 dxi dxi 



2 2 

Edu 3 duj du 3 s—^ du 3 du 3 du 3 
f)T ■ FIti FIti -<- — ' f)*r^ f)*ri finri 



E 



/ = ! 



du 3 
dx 3 

du 3 
dxi 



dui 
dxi 



3,1=1 

' du 2 
.dx 2 



dxj dxi dx, 



l OXi 



1=1 



dx 3 dxi dxi 



dx\dx2dx 3 



du\ du2 
dx\ dx2 



dui 
dx 2 



E 

.k=l 



du k du k 
dx 3 dxi 



E 

k=l 



du 3 du k du 3 du 3 
dx k dxi dx 3 dxi 



dx\dx2dx 3 



dui du2 
dx2 dx\ 



<C 



\u 3 \ |Vu| |V/jVm| dx\dx2dx 3 



Next, we will estimate the right-hand side of the above inequality using either Lemma 2 or Lemma 3. Each 
will be used for dealing with either one of the conditions (15) or (16). On the one hand by applying (13), with 



= M,/= \Vu\,g= \V h Vu\, and r = 



3a-2 



we get 



|m 3 | \Vu\ |V/jVu| dx\dx2dx 3 



<c\\u 3 \\P* \\d x M\l^\\^Af* 



R 2 v w ||r- 2 \\d x ^u\\i^ ||vv fcU || 2 



<c|| U3 ||2||^ lM3 ||°- i ||vu||°- 1 \\d x ^u\\^ + - || vv fc t 



(19) 
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On the other hand by applying (14), with <fi = |us|,/ = |Vu|,<7 = \\7 h \7u\, we obtain 
\us\ \Vu\ |VhVu| dx\dx 2 dx 3 



<C\\u^\Q 1 ^ \\d x ^\l^\\^u\\^ \\d Xl Vu\\f- 2 \\d X2 Vu\\f- 2 ||VV,,,7 



hU\\2 



<C|| U3 || 2 °- 2 \\d X3 uz\\Z- 2 \\Vu\\i + - \\VV h u\\i. 



In case we use (19) we obtain 
d\\V h u\{i 



<H 



^||V^Vm||1 < CllualllH^^all^- 1 HVwH^ ll^VtiHI"- 1 ; 



Alteratively, if we use (20), we obtain 



d\\V h u\\ 
dt 



+ v\\V h Vu\\i < C\\u 3 \\ 2 - 2 ||^3 M3 ||r 2 ||V u ||^. 



(20) 
(21) 
(22) 



Therefore, integrating (21) and using Holder inequality and applying (17) we get 
\\V h u(t)\\l + V ( \\V h Vu{ S )\\%ds 



ii 

j-t ..., \ 277 r 

|2 i r< I I IIP „. /„M|H=2 IIT7»./'oM|2 



t \ 2(c-l) 

2 



<llv h «oiis + c(/ II^ iU3 (s)||s- 2 i|Vm(s)||^ s ] u i|A«(«)i^d« 

for all * £ [0,7"*). Alteratively, integrating (22) and using Holder inequality and applying (17) we get a different 
estimate 

7u{s)\\tds< \\V h u \\l + cj (\\d X3 u 3 (s)\\f"Wu(s)\\i) ds 

for all* e [0,7"*). 

2.2. || Vu||2 estimates. Taking the inner product of the equation (1) with —Am in L 2 , we obtain 

2 



(23) 



(24) 



ld||Vw||j 
2 dt 



v\\ku\\ 



I [(u ■ V)w] • Afrit dx\dx 2 dx 3 + I [(u ■ V)u] • d X3 u dx\dx 2 dx 3 

JR3 JR3 



<c 



\u 3 \ |Vu| \V h Vu\ + \V h u\ \d X3 u\ x 



dx\dx2dx 3 



(25) 



(26) 



Applying the Cauchy-Schwarz inequality and (9), with r = 4, wc obtain 

|Vfc«| |dx 3 u| dxidx 2 dx 3 < ||Vh«||2 ||Vu||4 

< ||V fc «|| 2 ||Vu||£ /2 ||V fc V«||2 ||A U ||^ /2 . 

Now, we are ready to complete our proof for the case when condition (15) holds. By (19) and Young's inequality, 
we have 

\u 3 \ |Vu| |V/jVu| dx\dx 2 dx 3 
^CllwsllillSx^sll^llVullf 51 \\d X3 Vu\\f T + V - ||VV h «||l 



<C\\u 3 \\ 2 - 2 \\d xl u 3 \\i- 2 \Wu\\l + ^-\\Au\\l 



4(o-l) 
X-2 
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As a result of the above and (25), we get 

glgHJ + g ||An|£ 
dt 2 " " 2 



4(a-l) 
x-2 



< C||V h «|| 2 ||V M ||2 /z ||V ft V u || 2 ||A U ||^ + C||u 3 || 2 °- 2 ||^ lU3 ||«- 2 ||V u ||^ 
Integrating the above inequality and using Holder inequality we obtain 



\Vu{t)\\l+ V -j\\/±u\\lds 



<||Vu(0)||; + C( sup ||Vfc«(*)|| 2 

yO<s<i 



|Vw|| 2 ds 



|VfcV«(*)||| ds 



\\Au(s)\\ 2 2 ds 



C sup \\u(s) 



4(o-l) 
I a-2 



||^3 U3 ( S )||a- 2 ||Vu( S )||^ S 



Thanks to (17) and (23), we get 



l|V«(t) 



2 , ^ 



2 Jo 



Au|| 2 ds 



< || Vw(0) 1^ + Ci^ 



1/4 



VfcUolla + C / \\d Xl u 3 (s)\\£- 2 \\\7u(s)\\ids 



, ft 2q \ 2( Q -1) 

|2 , /-< / / II o „. /„\||H=2 I|r7„./„M|2 



\\&u{s)\\lds 



-CK 1 



\\d Xl u 3 (s)\\S- 2 \\Vu(s)\\Us 



By Young's and Holder inequalities, we obtain 

rt 



l|V«(t)||^ + 7 



o 



|Au|| 2 ds 



< C\\\7u(0)\\i + C (J \\d xl u 3 (s)\\S- 3 \\Vu(a)\\i ds 
"" ' ' ||9 XlU3 (s)||r 5 ||Vw(s)||^s 





\\Vu{s)\\lds 



2(c-2) 



Thanks again to (17), we get 



||Vu(£)|| 2 + - / ||A U || 2 ds<C\\Vu{0)\\ 2 2 + C \\d Xl u 3 (s) \\r 3 1| Vu(s)\\ 2 2 ds 
4 Jo Jo 

Therefore, in case (11) holds, we apply Gronwall inequality to obtain 



(27) 



(28) 



||V«(*)||2 + v / ||Aw|| 2 ds < C(l + ||V«(0)||2)e CM , 
Jo 

for all £ E [0,7"*). Therefore, if the condition (15) holds the H 1 norm of the solution u is bounded, and this 
completes our proof in this case. Next, we complete the proof when u 3 satisfies (16). Thanks to (25), (20) and 
(26), we get 

rfJIVwJj! v 2 

^^ + 2 I|AU| ' 2 

< C||V h u|| 2 ||Vn||^ /2 ||V ft Vn|| 2 ||Au||£ /2 + C\\u 3 \^ \\d x ,u 3 \\f~ 2 \F uf* 
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Integrating the above inequality and using Holder inequality we obtain 



\Vu(t)f 2 + ^J*\\Au\\Us 



<||V«(0)||; + C sup ||Vfc«(«)|| 2 



|Vu|| 2 ds 



\V h Vu(s)\\ 2 ds 



\\Au(s)\\l ds 



+C[ sup ||u(s)||j 



Thanks to (17) and (24), we get 



\\d X3 u 3 (s)\\£- 2 \\\7u( S )\\ 2 2 d S 



l|V«(t)||^ + ^ 



|Au|| 2 ds 



<\\Vum\l + CKll A 



Vhu \\ 2 2 + C( / \\d Xs u 3 (s)\\S- 2 \\Vu(s)\\Us 



|Am(s)||| ds 



+CK 1 



\\d Xl u 3 (s)\\r 2 \\Vu(s)\\ 2 2 ds 







By Young's and Holder inequalities, we obtain 



2 , y_ 

4 7o 



\Vu{t)\\ 2 + - I ||Ati||' ds 



+CK 1 - 
Thanks again to (17), we get 
l|V«(*)|£ + 7 



< c||v«(o)||^ + c (j n^^^lir- 2 ' ||v«(*)||^ ds 

^' ' ' \d X3 u 3 (s)\\f^\\\7u(s)\\ 2 2 ds 



|Vu(s)||| ds 



Au\\ z 2 ds<C\\Vu(0)\\i + C[ / \\d X3 u 3 (s)U la - 2) \\Vu(s)\\ 2 2 ds 



Therefore, by Gronwall inequality and (16) we obtain 



(29) 



(30) 



\\Vu(t)\\ 2 2 + v \\Au\\ 2 2 ds<C(l + \\Vu(0)\\ 2 2 )e CM . 
Jo 

for all t G [0, 7"*). Therefore, the H norm of the strong solution u is bounded on the maximal interval of existence 
[0, T*). This completes the proof of Theorem 1. 
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